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The universe displays a three-dimensional pattern of hot and cold spots in the radiation remnant 
from the big bang. The global geometry of the universe can be revealed in the spatial distribution 
of these spots. In a topologically compact universe, distinctive patterns are especially prominent 
in spatial correlations of the radiation temperature. Whereas these patterns are usually washed 
out in statistical averages, we propose a scheme which uses the universe's spots to observe global 
geometry in a manner analogous to the use of multiple images of a gravitationally lensed quasar 
to study the geometry of the lens. To demonstrate how the geometry of space forms patterns, 
we develop a simple real-space approximation to estimate temperature correlations for any set of 
cosmological parameters and any global geometry. We present correlated spheres which clearly show 
topological pattern formation for compact flat universes as well as for the compact negatively curved 
space introduced by Weeks and another discovered by Best. These examples illustrate how future 
satellite-based observations of the microwave background can determine the full geometry of the 
universe. 
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^From the zebra's stripes to the leopard's spots, the 
animal kingdom displays a diversity of coat patterns. 
Following the innovative ideas of Turing jlj , mathemati- 
cal biologists have posed and partly answered the ques- 
tion of how the leopard got its spots. The fluctuation 
of enzymes diffusing through the developing embryo can 
lead to the spatial pattern formation displayed by animal 
coats. Both the geometry and size of the animal exert a 
strong influence on differentiating patterns. For instance, 
the broad cylindrical shape of the leopard's body favors 
spots while the tapered tail induces stripes [|| . Remark- 
ably, these diverse features can arise from the proper- 
ties of simple solutions to second-order partial differen- 
tial equations on the geometry and topology appropriate 
for animal limbs, tails, or bodies. [] 

Similarly, the global geometry of the universe can lead 
to distinctive pattern formation by the normal modes of 
vibration of the universe in the microwave sky, or even 
in the distribution of luminous galaxies. Unobscured by 
complicated evolutionary effects, the cosmic microwave 
background radiation (CMB) offers the best site to seek 
out patterning. When light last scattered off hot matter, 
small temperature fluctuations left birthmarks in the ra- 
diation. These hot and cold spots in the primordial ra- 
diation may be randomly distributed in an infinite uni- 
verse, but if the universe possesses a compact topology 
then distinctive patterns can be generated. Such topolog- 
ically induced pattern formation has already been seen 
in simulated maps of the microwave sky 

The diffusive processes in biological systems actually 
simplify an animal's coat by singling out a particular 
mode and hence a particular pattern. By contrast, the 
universe contains many competing modes. While some 
simple large scale patterns are clearly evident in sky maps 
of the CMB temperature fluctuations (|||, subtler pat- 
terns can be extracted from the observations. In par- 
ticular, a correlation in temperature between pairs of 
points can scan the map and select out geometric fea- 
tures. Maps of correlations for different topologies pop- 
ulate the cosmic zoo of possibilities displayed in sections 
0andg 

CMB observations of the universe's hot and cold spots 
could be used as real-space pictures of geometry. With 



correlated maps we could observe topologically lensed im- 
ages of the horizon at the time of last scattering. The 
number and pattern of lensed images of the horizon pro- 
vide information about the geometry of the universe. A 
nice analogy is provided by observations of gravitation- 
ally lensed images of a quasar, the pattern and num- 
ber of which reveal the geometry of the intervening lens. 
The universe's spots are too small in angular scale to 
have been detected by the Cosmic Background Explorer 
(COBE) but will be visible by the future satellite mis- 
sions, MAP (Microwave Anisotropy Probe) and Planck 
Surveyor. 

The fluctuations in the CMB can be described 
by decomposing the relative temperature fluctuations, 
ST(x)/T — J d?k(j)pjj^ into a series of eigenmodes ip^ 
each of which is a solution of Laplace's equation on the 
curved space of the universe, (V 2 + k 2 ) = 0. The 

are primordially seeded amplitudes that on average 
define a spectrum of perturbations. Topology intro- 
duces boundary conditions which create a discretized set 
of wave vectors and complex relations between the (j)^ 
[|[. This problem is tractable in the case of a flat uni- 
verse but, for a compact hyperbolic space, the identifica- 
tion rules lead to such intricate boundary conditions on 
the modes that they cannot be decomposed analytically 
[f?J. This is a symptom of the chaotic hyperbolic flow 
of geodesies induced on these spaces. Microwave pho- 
tons moving on compact hyperbolic manifolds thus de- 
mand more indirect methods of analysis. Cosmological 
observations favor a negative curvature and mathemat- 
ics favors compact hyperbolic manifolds by providing an 
infinite number of candidates. The resistance encoun- 
tered when subjecting these manifolds to conventional 
eigenmode treatments has sparked recent interest and 
produced several different approaches to the search for 
observational probes of the universe's global topology. 
Three strategies have emerged: a search for circles in the 
microwave sky [|j , a search for spatial pattern formation 
[UH, an d a direct attack via the method of images [||. 
All three approaches are necessarily related. We will see 
aspects of each emerge in our zoo of topological examples. 

For reviews and other interesting papers on cosmic 
topology see Refs. [foj-OI. 



1 For example, the two-dimensional Helmholtz equation 

V 2 + /c 2 = O, 

with (n • V)0 = on the boundary, has solutions which 
describe different geometric tesselations of alternating re- 
gions with <j> > and (f) < 0, For example, an infinite 
line of one-dimensional stripes is described by the solutions 
<j) — cos kx, with k = nir for n — ±1, ±2, .., whereas a solution 
cj) = \ [cos kx + cos ky] , with k = tv or 2ty describes a chequer- 
board tesselation of space with alternating square spots, with 
cj) > and <j> < 0, inclined at 45 o . 



I. TOPOLOGY AND THE ANGULAR 
CORRELATION FUNCTION 

Compact spaces are specified by Ai = U/T where U is 
the geometry of the infinite space before any compactifi- 
cations are made. U is also known as the universal cover- 
ing space. The generators of the group Y provide the 
set of instructions for gluing together points in U in or- 
der to render the space finite and multiconnected. The V 
forms a discrete subgroup of the full isometry group of the 
covering space. A finite geometry has the same curvature 
as U. The universal covering space is usually assumed to 
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be a space of constant curvature, either the flat E 3 , the 
negatively curved H 3 , or the positively curved S 3 . The 
real universe does not have perfectly constant curvature 
over the entire manifold: its true topologies will then be 
a deformation of perfect polyhedra. The statement that 
the universe looks homogeneous and isotropic on average 
indicates that constant curvature will be a reasonable 
first assumption for finite manifolds which are compa- 
rable to or smaller than the Hubble scale of the visible 
universe. We limit ourselves to the subset of observation- 
ally viable manifolds with covering spaces of either the 
flat E 3 or hyperbolic H 3 variety. 

In order to scan the sky for evidence of compact topol- 
ogy we use the correlation function for the temperature 
at two different points on the sky, 



ST(n) ST(n f ) 



= C(n. n ). 



(1) 



where ft is a unit vector pointing from the Earth towards 
some location on the sky. The sphere of radius Ar]h de- 
fines the surface from which we receive last-scattered pho- 
tons. The size of the radius, A77, is the conformal time be- 
tween the time of last scattering and the present. The an- 
gular brackets denote an average over all possible realiza- 
tions. In a simply connected cosmos, the assumption that 
the hot and cold spots are homogeneously and isotropi- 
cally distributed ensures that the correlation function de- 
pends on only one parameter, namely the angle between 
the two points on the sky, that is C(n, n') = C(0) where 
cos# = n • h' '. All of the information in an homogeneous 
and isotropic Gaussian field is contained in C(6). 

By contrast, topological identifications always break 
isotropy (with the exception of the projective space on 
S 3 ) and most break homogeneity as well. The hypertorus 
on E 3 is the only homogeneous compact space. Conse- 
quently, the correlation function, C(n, 77/), depends fully 
on n and n' and would require four dimensions for a full 
representation. Predicting C(n, n') is also more challeng- 
ing on compact manifolds, particularly the most unusual 
compact hyperbolic spaces. 

We will introduce a real-space approximation which 
captures all the important features necessary for examin- 
ing global structure. Our approximation is based on two 
physical observations. Firstly, two points which appear 
to be widely separated may actually be close together, as 
in Fig. |l|: there is not a unique distance between points. 
Secondly, the correlation between two points is strongly 
peaked around short separations. While global topology 
drastically changes the large scale perturbations, it does 
not have a strong impact on scales much smaller than the 
size of the physical space. Therefore we expect the corre- 
lation function between nearby points to be well approx- 
imated by the correlation function of a simply connected 
universe. 

Motivated by these two observations, we propose the 
following real space approximation to the angular corre- 
lation function: we take the correlation in temperature 
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FIG. 1. Two points may appear to be far apart but if we 
identify opposite sides of the rectangle we see that the same 
points are actually close together on the compact torus. 



between two points on the surface of last scattering to be 
the correlation function in a simply connected universe 
given their minimum separation, that is, 



C M (n,n) « C u [d min (x(h),x(n))} 



(2) 



where C u is the correlation function on the universal 
cover, x(h) is the physical location of a point on the sur- 
face of last scattering, and d m [ n is the minimum distance 
between the two points in the topological space. To find 
the minimum distance we locate the image points with 
the generators of the identifications. The point x' has 
first neighbor images at locations = g^x' and second 
neighbor images at locations y^kx — gk 2 9k 1 %' \ etc. More 
concisely, the images out to order m can be written as 



X\gk,x'{ti) 



(3) 



The image point which lands closest to x(h) determines 
^min- Note that all of the effects of topology appear only 
in the minimization of the separation between the two 
points. Thus, the correlation function can be approxi- 
mated using only a knowledge of the identifications that 
compactify the space: the very relations that are used to 
specify the topology. 

Our approximation is closely related to the method of 
images employed in Refs. The correlation function 
computed by the method of images is a sum over all of 
the copies of the image points in the simply connected 
space with the same curvature, 

C(n,n') = lim V C u {d(n, 7*n')) 

r* ^oo — ' 

li 

"17 I""* drsmh 2 rC u (r) , (4) 
V Jo 

where the second term is a regularizer introduced to con- 
trol the exponential proliferation of images, particularly 
on H 3 , as the radius r* increases relative to the volume 
of the fundamental domain, V The 7 are composite 
elements of the group T. Our approximation amounts to 
keeping only the dominant term in the sum over images. 
Maintaining the first term is by far less cumbersome than 
summing an infinite number of terms and fares well in 
approximating the exact correlation function. The ap- 
proximation is also quite valuable for our purposes since 
we are able to easily include the essential effects of the 
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smoothing across the horizon at the time of last scat- 
tering and the relevant microphysics at work over small 
separations. The small-scale physics at the time of last 
scattering is important for distinguishing topologically 
l ensed images from fictitious correlations as explained in 
B[ The physical processes operating at last scattering 
can in principle be folded into the method of images to 
refine our approximation and sharpen the focus on our 
pictures of geometry. The thickness and velocity of the 
surface of last scatter will also induce additional correc- 
tions. A test of our approximation when applied to the 
flat geometries is given in detail in § |II B , 

Given an estimate of the full four-dimensional 
C^^n'), we could build realizations of a map of 
5T(n)/T for any topology, even without the eigenmodes. 
A given universe would be obtained by a random real- 
ization of a Gaussian-distributed variable with a mean 
(ST(n)/T) = and variance Cm(ti, n'). However, we 
can do better. As Cm{^^') is a function on two copies 
of the sky, it necessarily must be calculated numerically 
for Np ix values, where N pix is the number of pixels into 
which the sky is divided. This is cumbersome, if not 
practically impossible, to carry out for many topologies 
on a reasonably finely grided sky map, particularly if we 
aim at the resolutions expected of future CMB satellite 
missions. Furthermore, many of the peaks of this func- 
tion occur at values where ft and ft' are separated by a 
small angle on the sky. Since these points would be close 
together regardless of global topology, much of the base 
space of Cm{^^') is useless in discriminating between 
topological properties of cosmological models. 

Instead of building full ST(n)/T maps, we build se- 
lective correlation maps. An illustrative example is the 
antipodal map constructed by evaluating the correlation 
function at antipodal points along the sky, 

k{n) = C M {fi,-n). (5) 

In a simply-connected universe A(n) would produce on 
average nothing more than an overall monopole. Fortu- 
itous correlations may appear at random in a given re- 
alization but there would be no defining structure. In a 
universe with compact topology a great deal of structure 
can be surveyed by means of this one simple correlation. 
While the best analysis for any particular topology will 
make use of the full correlation function, the antipodal 
map will prove to be a most useful survey tool. To get a 
feel for the way in which correlations reflect global topol- 
ogy, we make a gallery of correlated spheres first for all 
flat topologies and then for a few choice compact hyper- 
bolic spaces. We also consider correlations under symme- 
tries other than antipodal as the specific space demands. 
In order to represent the correlations as 2D maps we 
consider correlations of the form 

C g (h) = C M (n,gn/\\gh\\) , (6) 

where again g is an element of the group T. For the 
compact hyperbolic spaces, we also correlate an arbi- 
trary point with the rest of the surface of last scatter. 
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FIG. 2. The hexagonal geometry. The prism face is glued 
with a twist of 2tv/3 to create a topology distinct from the 
hypertorus. 

In many ways these point-to- sphere correlations, are the 
most promising. They very dramatically reveal geomet- 
ric patterns and they do not require any foreknowledge 
of the symmetries of the space. 

II. THE COSMIC ZOO 
A. Flat Topologies 

Compact flat spaces have already fallen out of fa- 
vor as small universes. The first suspect for investiga- 
tion was the simplest hypertorus Q. The fundamen- 
tal domain is a parallelepiped with opposite faces iden- 
tified in pairs. As shown by Stevens, Scott and Silk 
Jl5[ , the square hypertorus suffers a sharp truncation of 
long wavelength power in temperature fluctuations, too 
sharp to be consistent with COBE observations unless 
the box is larger than about 40% of the observable uni- 
verse. Later, anisotropic tori were studied in Refs. Q 
using symmetry methods, and were similarly bound. The 
tightest limits on an equal-sided square hypertorus were 
obtained in Ref. (|] using the method of images, placing 
the topology scale just beyond the observable universe 
with h > 2.19A?7 where h is the length of the side of the 
square and A77 is the conformal time since last scattering. 

With the universal covering space of E 3 , there are only 
five more orient able, compact spaces that can be con- 
structed ]l8[ . Three are built by identifyin g th e faces of 
a parallelepiped with relative twist (see Fig. |l0|), and two 
are constructed by identifying the faces of a hexagonal 
prism with the prism faces twisted by 2tt/3 or tt/3 relative 
to each other before being identified. In Ref. we 
derived the eigenmode spectrum explicitly for all twisted 
cases and showed that the angle- averaged power spec- 
trum is incompatible with COBE for equilateral spaces 
unless the universe was really big, if not actually infi- 
nite (i.e. topological identification scale > 80%A?7). The 
compact flat spaces are thus of limited interest, although 
polyhedra with sides of disparate lengths may still be 
viable. 

However, since we do have eigenmodes for all flat 
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FIG. 3. A typical map of 8T/T(fi) in a hexagonal prism 
with h O.8A77. 



topologies |J, these topologies are important testing 
grounds for any method which attempts to circumvent 
the eigenvalue problem. We therefore use E 3 to test our 
real space approximation before moving on to compact 
topologies of H 3 . 

For illustration, we spend time explaining the features 
in the markings for a 27r/3-twisted hexagonal prism. The 
manifold is E 3 /r with T the group of instructions for 
identifying the faces of the hexagon. Define a coordinate 
system for which e\ is orthogonal to the hexagonal face 
and the non-orthogonal vectors £2, £3 span the face as in 
Fig. [| T is generated by #2? #3} where g\ generates a 
rotation through 2n/3 about the e\ axis combined with 
a translation orthogonal to the hexagonal face through 
a distance c, #2 effects a translation along 62 through 
h, and gs effects a translation along is also through h. 
Another way to visualize E 3 /r is to glue copies of the 
fundamental domain together according to the identifi- 
cation rules. In this way E 3 can be completely tiled with 
layers of hexagons separated by the length of the prism 
direction. 

We can use these symmetries to build the discrete spec- 
trum of eigenmodes which describe fluctuations on this 
space as was done in Ref. Using these modes a typ- 
ical map of ST(h)/T for a simulated compact hexagonal 
universe can be created. We show this for comparison 
in Fig. ||. There is something at work in the map of 
ST(n)/T but it is hard to define. 

If, however, we inspect the idealized zero- variance an- 
tipodal map, the pattern jumps out at us. These antipo- 
dal maps were built with our approximation and did not 
require the full mode solutions. We locate a point on the 
sphere of last scattering (SLS) x — Ann and compare 
the correlation of that temperature with its antipodal 
point x' = —Ann'. If the universe were simply E 3 , that 
is, flat and unconnected, then these points should be to- 
tally uncorrelated on average. Given that this space is 
multiconnected, the two opposed points may in fact be 
close together as demonstrated for the 2D surface of Fig. 
p]. We approximate the correlation between antipodal 
points as the correlation that the points would have in a 
simply connected space given their minimum separation 
as in (0). Antipody is then approximated as 




FIG. 4. C u (0) for a flat, COBE normalized CDM cosmol- 
ogy. Note that the becomes negative due to the definition of 
ST/T. 




(A(n))~C u [d min (x(n),x\-n))] 



(7) 



FIG. 5. Orthographic projection of A(n) at a resolution of 
20 arcminutes for a hexagonal prism with h — b — c — O.8A77. 
The observer is at the origin. 



where C u (0) is obtained from CMBFAST @ for a stan- 
dard, flat CDM cosmology (Q b = .05, ^ c = .95, H = 
SOKms -1 Mpc -1 ) with the dipole component calculated 
as per a flat power spectrum. This is illustrated in Fig. 
||. We minimize the distance by first taking the image 
of these points under the action of the generators of the 
group T until we relocate them within the fundamental 
domain. Leaving one of these relocated points fixed, we 
consider all images of the second point that lie within one 
of the nearest neighbors of the fundamental domain, and 
choose d m i n as the shortest distance from the first point 
to one of these images. Note that only by considering 
all of these images, including those that are diagonally 
located relative to the fundamental domain, can we be 
sure that our definition of distance depends on the overall 
topology of the space rather than on the particular coor- 
dinates that we use to fix the fundamental domain. This 
method results in an overall monopole component of the 
antipodal map which we wish to discard. We simply re- 
move the monopole as calculated from the antipodal map 
alone, with the understanding that this is to be compared 
with a measured map that is similarly normalized. 
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FIG. 6. The surface of last scatter intersects the layered 
tiling of flat space. Each full tile represents a copy of the 
fundamental domain. Slice 1 is represented on the lower left 
and slice 2 on the lower right. The dots show correlated points 
picked up in the antipodal map and explain the emergence of 
the hexagonal geometry in Fig. ^ 

It is customary to use the Aitoff projection (as in Fig. 
||) to view the map of ST(h)/T so as to see the entire 
sky. For the antipodal map we prefer the orthographic 
projection which shows the genuine shape of the surface 
of last scattering. In this three-dimensional view, there 
is no distortion of the sky pattern. No information is lost 
by limiting ourselves to a view of only half the sky at a 
time since A(n) is by definition symmetric under n. 

There is clearly a hexagon in the antipodal map. The 
intersection of the hexagonal layers with the spherical 
surface of last scattering defines rings centered on the 
prism direction and are clearly picked out by A(n). The 
rings of correlated spots occur since the correlation be- 
tween points on the sphere of last scattering separated 
by 7r for this geometry is the same as the correlation 
between opposites sides on a circular slice through the 
sphere taken along the hexagonal plane, as demonstrated 
in Fig. |[ The bright primary spots are identical points 
at the core and are just correlated away from the center. 
In addition to these identical images, there are also sec- 
ondary spots picked up due to the correlation of regions 
which just near each other. 

These rings of structure are not the circles of the sky 
of Ref. (|], although they are related. As pointed out 
there, pairs of identical circles occur in the microwave 
sky due to the intersection of the surface of last scatter- 




FIG. 7. Orthographic projection of A(n) for a hexagonal 
prism with h — b — c — O.6A77. 

ing with copies of itself. None of the circles are located 
by A(n) for the 27r/3 twist because none of the circles 
in this space are paired under a 7r-symmetry. For rings 
separated along the prism direction by multiples of 3c, 
the hexagonal faces have completed a full rotation and 
antipody compares one point to its opposite face. We 
continue to call the concentric collections of spots in the 
correlated A{n) maps 'rings' and we reserve 'circles' for 
intersections of the copies of the surface of last scatter- 
ing with itself. So, the primary spots on a given ring 
lie on different pairs of circles in the sky. The rings of 
secondary spots do not lie on circle pairs. 

As the topology scale gets smaller, there are more rings 
of patterns and the first ring appears closer to the caps 
as shown in Fig. [?]. 

The size of a spot will be set by the Silk damping 
which smooths fluctuations on small scales. For sepa- 
rations which exceed this length, we expect correlations 
quickly to die off. Since the Silk damping scale is smaller 
than the horizon size at the time of decoupling we expect 
the angular size of these spots to be too small for COBE 
to have detected. The beam smearing in the COBE ex- 
periment would dilute a spot over such a large angle that 
these bright markings would be below the signal to noise 
sensitivity of the detector. However, any experiment 
with high enough resolution and sensitivity to probe the 
Doppler peaks will be able to resolve these pattern on 
the sky. The planned missions MAP and Planck Sur- 
veyor would thus be able to detect the universe's spots. 

The other topology built from the hexagonal plane 
tiling involves a relative twist of the prism face through 
7r/3. As seen in the left-most panel of Fig. ||, there are 
again correlated spots on rings. When we make the space 
sufficiently small so that at least three copies of the fun- 
damental domain fit within the SLS the first pair of cir- 
cles appears in A(n), as shown in the right-most panel of 
Fig. I 

The other four topologies are built from a paral- 
lelepiped as the fundamental domain, as for the simplest 
hypertorus. No circles will be located by A(h) for the 
hypertorus, which is not to say that the hypertorus has 
no circles, just that none of the pairs of circles are ro- 
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FIG. 8. A(h) for a hexagonal prism with a 7r/3 twist 
h — c — 0.75A?7 there are no circles. To the right the z 
direction is .24 while h = 1. There are circles 




FIG. 9. A(h) for the torus with h=b=c= O.3IA77. 

tated by 7r relative to each other. Since the hypertorus 
does not involve any rotations of the faces, all three di- 
rections are picked out equally by antipodal pairings, as 
seen in the small torus of Fig. ^. The square geometry of 
the fundamental domain is blatantly traced out by the 
correlated spots leading to the disco ball effect. 

For the twisted parallelepipeds, A(n) will locate the 
axis of symmetry along which the faces are twisted. The 
correlated spots still trace out the symmetric square of 
the equilateral untwisted directions, as shown in Fig. [h]. 
For the ir/2 twisted space, the circles appear in A(n) 
if more than two copies of the fundamental domain fit 
inside the observable universe. At least one pair of circles 
will appear always for the 7r-twisted space as well as for 
the last compact topology built by gluing a parallelepiped 
by a series of diagonal translations and tt twists. 

The menagerie in this flat zoo is a testament to the in- 
fluence of geometric patterns on the universe's markings. 
An antipodal correlation is just one statistic but note 
how clearly antipodal correlations reveal the other sym- 
metries of the space. For instance, the prism direction 
perpendicular to the hexagonal face is clearly identified. 
Once the symmetries of the space do begin to become 
apparent, the future data on the CMB can be systemati- 
cally scanned for correlations under other symmetries of 
a given topological space: Cm (n, gh/\ \gh\ |). In this way 
the entire fundamental group T can be isolated, thereby 
determining the full geometry of the space, as was also 




FIG. 10. A(h) for a 7r/2-twisted space. 




FIG. 11. The sphere of 

correlations C(h(x,y, z),h(x,y, —z)) for the hypertorus. The 
topology scale is O.4A77. 

suggested in Ref. |2(| . 

The indicator A(n) is particularly convenient because 
it can be visualized on an half-sphere without the specifi- 
cation of a particular axis, but once an axis of symmetry 
is chosen, correlations at other angular separations can 
similarly be examined. As an example of such a corre- 
lation, we consider Cm{^^') where we merely compare 
the point (x, y, z) to the point (x, —z) for the hyper- 
torus of size 0.4 of the SLS. We assume that this axis has 
been identified by examining the overall antipodal prop- 
erties of the sky. We expect this correlation to pick up 
circles and it does, as demonstrated in the Fig. [ll[ 

As another example we compare, for the 7r/2-twisted 
space of size 0.4 the radius of the SLS, the correlation 
between two points separated by a rotation of tt/2 around 
the x axis. In Fig. [L2] we display the result of applying 
the correlation function 

C M [drnin{x{0,4>),x'{0^ + K/2))] . (8) 



B. Observing the Universe's spots 

Our topological zoo of maps makes a compelling visual 
argument for a pattern-driven approach to the search 
for cosmic topology. Nevertheless, they are still only 



FIG. 12. The correlated sphere comparing points on the 
surface of last scatter with points related by a 7r/2- twist 
around the z-axis. The space is a 7r/2-twisted square with 
sides of length O.4A77. 



ensemble- averaged images calculated using an approxi- 
mation method. As the real microwave sky contains no 
such averages or approximations, we are left wondering 
if we should believe our eyes. We must look more closely 
at what is observed. 

In this section, we address the accuracy of our method 
by comparing our approximate maps with exactly calcu- 
lated ensemble- averaged maps as well as with simulated 
realizations of the microwave sky. We consider the fiat 
hypertorus which has a simple eigenmode decomposition 
of the temperature fluctuations, 
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(n) oc 0£ exp (^iAr]k • fi^j , 



FIG. 13. A comparison of the real-space approximated an- 
tipodal map (left) with the exact ensemb led- averaged map 
(right) for a torus space with h — b — c — .8A77. 




(9) 



FIG. 14. Real-space approximated antipodal map, cal- 
culated from a flat power-spectrum in which the dipole, 
quadrupole, and octopole have been removed. 



with k = 2ir(n x /h,n y /b,n z /c). The (j)^ are primordially 
seeded Gaussian amplitudes that obey the reality condi- 
tion = </>* ^ and as an ensemble define the spectrum 
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(10) 



With this decomposition we can construct the correlation 
function between any two points on the sky as 



C(n,n') = ('^(n)^(n') 



(11) 



oc 



k 3 



■ exp 



(iArjk • (n — n')J . 



The antipodal correlation on the hypertorus is the simple 

case, C(n, — h), and 



(A(n)> oc V exp(i2Arjk ■ 



n). 



(12) 



up to an overall normalization. This is nothing more than 
a Fourier transform of the power spectrum, and can be 
computed exactly. 

In figure [l3| we compare the results of our approxima- 
tion with the exact eigenmode decomposition of eqn. (O) 



for the case of the hypertorus with h = b = c = 0.8. In 
both calculations, we take a simple flat (V(k) = 1) power 
spectrum with no additional physics added. 

In both maps the small-scale structure and location 
of the spots are identical, but in the approximate map 
the large-scale features of the spots are smeared out. 
This smearing stems from the fact that our approxima- 
tion does not incorporate the damping of power at low 
wavenumbers due to the discretization of eigenmodes on 
the compact space. Thus C c/ [<i m i n (x(n), x'(h'))} falls off 
more slowly than the corresponding correlation function 
on the compact space. 

This discrepancy can be remedied by replacing the 
C u with an angular correlation function from which the 
lowest multipoles have been removed. In the case of a 
h = b = c = 0.8 space, the damping is well approximated 
by removing the dipole, quadrupole, and octopole terms 
[|l7f . This results in the approximate map shown in Fig. 
O , which is quite close to the exact result. In principle, 
the broadening of the spots can also be removed by in- 
cluding more terms in the method of images expansion, 
eqn. ||. In practice, however, a cumbersome number of 
terms have to be included in order to begin to approach 
the exact result, with the maps actually getting worse 
before they get better. 

The question remains how to analyze real data so as 
to extract patterns with confidence. It is customary in 
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an analysis of CMB observations to take angular aver- 
ages in Fourier space. All standard methods thus smear 
out the very patterns we seek. Instead we advocate that 
the correlation maps be treated as real space pictures of 
geometry. As such, they are akin to pictures of a galaxy 
or to observations of gravitationally lensed images. 

All of the previous correlation maps are ensemble aver- 
ages. Since we have only one universe to observe and only 
one realization of the data we might worry that cosmic 
variance would drown out any of the features of topol- 
ogy in a given realization. We are unable to combat 
cosmic variance by averaging over the sky since it is pre- 
cisely such averaging which we are trying to avoid. We 
have numerically simulated high-resolution observations 
of a toroidal universe to demonstrate what future satel- 
lite observations will offer. Cosmic variance is not a ter- 
rible hindrance and we are able to extract the correlated 
information from the high-resolution simulations. 

The map of 5T(n)/T at the top of Fig. [IB] is a simula- 
tion of an all-sky map generated using the explicit eigen- 
modes of @. The physical processes at work on very 
small scales were modeled so as to emulate the Doppler 
peaks of a flat CDM universe. To obtain fairly accurate 
maps we have to keep account of the physical processes 
at the epoch of decoupling which determine the shape 
of the spectrum. At decoupling, the comoving length 
of the horizon is r/dec — 10 -2 A?7, which is presumably 
far smaller than the dimension of our fundamental do- 
main. The causal processes that modulate the spectrum 
are therefore not modified by the compact nature of the 
space. In this case we can use a suitable function for 
T(k) which is able to reproduce the height and position 
of the first Doppler peak in a flat CDM universe. The 
function we choose is the following: 



(^ n - 1 + n n+1 + ^ n+3 ) exp (-^ 2 /4) (13) 



where n = 10 -2 /cA?7 and n is the primordial fluctuation 
spectral index. We show this function in Fig. [l6[ In 
the same figure we also show the /c-range covered by 
the model we examined. With the power spectrum of 
eqn. (|l3|), the relative height of the first Doppler peak 
and the Sachs- Wolfe plateau are in good agreement with 
a flat CDM model for 0.8 < n < 1.2. As reported in 
H , the lower limit on k is fixed by the dimensions of the 
fundamental domain, while the upper limit is constrained 
by the resolution grid in k-space used to generate a given 
realization. We used a grid of 

Tlx f^y z 

= 256 3 elements in 

k-space which we fast-Fourier transformed back into real 
space to calculate our ST(n)/T spectra. The maximum 
wavenumber for a given realization is k max = n x ir/h. It is 
apparent from Fig. [L6]that it is unnecessary to go beyond 
this resolution limit to probe down to the accuracy of our 
approximate power spectrum. 

To obtain the angular power spectrum, Ct = 
^2 m \ a ^m\ 2 + 1) 7 we decompose the temperature 
fluctuations into spherical harmonics, 5T(h)/T = 




FIG. 15. The unsmoothed numerical CMB data with 
1024 x 513 pixels for a small square hypertorus with side of 
length O.2A77. The power spectrum for this simulation. 




100 1000 

kAr] 



FIG. 16. The power spectrum as a function of kAr]. 

X^m a irnXim{n) , using a modified version of the fast code 
developed by Muciaccia et al. J2^]. The power spec- 
trum of the simulation can be seen at the bottom of Fig. 
[l5| . Notice the enhancement of power due to the multi- 
ple copies of the fundamental domain, a feature already 
noted in Ref. fl7|] . While the C/s certainly do not con- 
tain all of the information in a map of ST/T in a universe 
with multiconnected topology, the spectrum does reveal 
the essential behaviour of 5T/T at very high £s where 
topology is less influential. 

The leftmost panel in Fig. [l?] is the antipodal map 
from the data smoothed on scales of 1.5°, as though the 
data were convolved with an experimental beam. The 
ensemble average computed from eqn. ( |l2|) is shown in 
the rightmost panel of Fig. [l7[ Clearly, the realization 
shows the structure of the ensemble average. While this 
small universe can be ruled out with the COBE data due 
to the conspicuous lack of power on small scales [p^|,p[, 
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FIG. 17. A comparison of the antipodal map read directly 
off the simulated data shown in Fig. ^ after smoothing on 
scales of 1.5° (left) with an ensemble average (right). 



the pattern would go undetected in the existing satellite 
data. If we had smoothed on scales of 7° to emulate 
the COBE experiment, any definitive pattern would have 
been washed away. The much higher resolution of MAP 
and Planck Surveyor is required to measure the universe's 
spots. 

The very small space was chosen for dramatic effect. 
Since these spaces sample low £ modes more sparsely, 
there is less cosmic variance in some sense and the pat- 
terns are easier to detect. The large spaces can suf- 
fer more contamination from low £ modes. To recover 
the bare information in the ensemble averages, the low £ 
modes may need to be cleaned from the larger spaces. 

We also consider the correlation which compares z — > 
— z, as this is a true symmetry of the space. If we align 
the fundamental domain with the z axis of the SLS, then 
the correlation compares points in the direction n(6, (j)) 
with points in the direction h'fjr — 0, <j)): 

c z (h) = c(h(o^),ti(7r-e^)). (14) 

Notice that even in a universe with no topology there 
will be some structure in such a correlation function. 
The equator so defined is always compared with itself 
and so shows more correlations than the poles even 
without multiconnected identifications. Since this is an 
actual symmetry of the hypertorus we should find 10 
circles in one hemisphere for the small space of size 
h = O.2IA77. We pick up nine circles at latitudes 
64°, 53°, 44°, 37°, 30°, 24°, 17°, 11° and 6°, and ten if we 
include the one identified point at the poles, as shown 
in Fig. [l8| (compare with Fig. [y]). It is difficult to see 
that these spots do in fact lie on fully correlated circles. 
After all, the maps only represent the simple product 
ST(n)ST(h f )/T 2 . A measure of the correlations across 
the pairs of circles would draw out the feature more 
crisply (see also, in this connection, the statistic sug- 
gested in Ref. ||). To illustrate, we plot the quantity 

Ufl)= ^ STjnie^)) 5T{n\K-e,4>)) ^ 2 ^ 
for the unsmoothed data, in the rightmost panel of Fig. 




FIG. 18. On the left is the correlated sphere C z (n) defined 
in eqn. [l4| read directly off the numerical data. On the right 
the quantity £z(n) of eqn. is plotted to draw out the circles. 




FIG. 19. On the left is C z (fi) read directly off the numerical 
data for the larger torus of size 0.41A?7 smoothed on scales of 
0.5° . On the right £ z (n) of eqn. is plotted. 



|l8| . This measure of the temperature difference singles 
out the circles. 

For a larger square torus of size h = O.4IA77, the cor- 
relation C z (n) of eqn. ( |l4| ) also finds faint circles. Fig. 
|l9| shows the occurrence of circles in the sky at latitudes 
of around 55°, 38°, 24° and 12°. The left figure is a map 
of the simple correlation function C z (n) of eqn. (|l4) 
smoothed at 0.5° while the right figure locates the circles 
more distinctly, without smoothing, by plotting the £ z (n) 
of eqn. (|l5|). The thin dark circles are the collection of 
identical points for which ^ z = 0. 

For contrast, we compare the predictions for a hyper- 
torus with those for a simply-connected flat cosmos. Fig. 
|20| shows a realization of the SLS for a flat CDM universe. 
The data is again smoothed on scales of 1.5° . The low £ 
modes missing from the map in Fig. |l5]are clearly present 
in an infinite universe. There is no evidence of a pattern 
in antipodal correlations, nor in the C z (n^n / (z — > —z)), 
as demonstrated in the bland pictures of Fig. [20|. 

In general we suggest a two-step investigation to ob- 
serve topological lensing as the data becomes available. 
Firstly, it is advantageous to smooth the data on some 
small scale. Correlations can then be read off more eas- 
ily from the smoothed map. Secondly, the unsmoothed 
data can then be exploited to measure the temperature 
variation across the spot at the location indicated by the 
correlated signal. The measured spectra of the related 
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FIG. 20. Simulated CDM data for a universe with no topol- 
ogy. The standard map of the surface of last scatter is shown 
on the top. On the bottom left is the map A(h) and on the 
bottom right is the map C z (fi). 



spots can then be compared. If the spectra match, topo- 
logically lensed images have been measured: a correlation 
measure across two spots can be used to confirm lensed 
copies of a fluctuation. The small-scale structure of the 
two regions should be almost identical,^], although large- 
scale effects such as the subtraction of the CMB dipole, 
or the integrated Sachs- Wolfe effect in open spaces may 
cause a difference in the mean temperature of the two 
patches. In practice, a generalization of eqn. ( |l5| ) such 
as 

^ = ^r(£-£ spo t) 5T(y-y 8pot ) y ^ 

where the temperature fluctuations are defined with re- 
spect to the local mean temperatures of the spots, should 
do well in confirming lensed copies of fluctuations. 

As we improve our understanding of the correlation 
function on short scales we might determine the spectrum 
of each of these spots and use this information to distin- 
guish fictitious correlations from the real thing. In an ac- 
tual realization, all of these correlated spots will be there 
under a web of spurious correlations due to cosmic vari- 
ance. These random correlations can be distinguished 
much as foreground sources are distinguished from the 
galaxies they occult. Spurious correlations will not be 
distributed on rings, nor will they have the characteristic 



2 The sphere of last scatter cuts through the same 3D volume 
differently, so even if a point is identical to another on the 
SLS, the surrounding small region can be a slightly different 
sample of the 3D patch. 



size and spectrum of the topological correlations. Cosmic 
variance as a form of cosmic noise could thus in principle 
be subtracted off the maps. Again, the task is similar in 
spirit to distinguishing the gravitationally lensed images 
of a quasar from other unrelated or foreground sources. 
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III. THE HYPERBOLIC ZOO 

We turn now to the application of our correlated 
spheres to compact hyperbolic universes. Compact hy- 
perbolic spaces are inherently chaotic. The exponentially 
deviating trajectories of geodesic motions on a space of 
negative curvature mix and fold chaotically through the 
space as they exit and enter the multi-faceted fundamen- 
tal domain. Chaos endows these spaces with many in- 
triguing properties, including fractal structures within 
the pattern of entangled geodesies. The patterns in- 
scribed in these skies thus promise to be intricate. Pri- 
mordial quantum fluctuations which ultimately seed the 
hot and cold spots are described by quantum chaos for 
which there are very different predictions than for non- 
chaotic quantum systems [^2|. The assumption of a flat, 
Gaussian seeded spectrum of perturbations may be a 
poor one. 

While complicated, these spaces are not obscure. 
There is a countable infinity of topologically distinct com- 
pact H 3 spaces, although they have yet to be completely 
classified |23|]. Furthermore, observations favor a universe 
with subcritical density. There is therefore considerable 
interest in understanding the predictions for these spaces. 
Regardless of how predictive maps are produced, the ul- 
timate question is, how do we analyze the CMB data to 
search for topology. Canonical treatments rely on angu- 
lar averages which smear out patterns and a likelihood 
analysis based on Gaussian statistics. While this may 
have some restricted meaning, it is dangerous to draw 
precise conclusions from a Gaussian isotropic probability 
distribution when the space itself destroys isotropy and 
the primordial spectrum is unknown. Furthermore, it re- 
quires a case by case analysis and may even depend on 
the location of the observer. 

This implies that a statistical analysis of the data re- 
quires a model template. If the universe is not a per- 
fect manifold of constant curvature the template match 
is lost. Instead of asking the statistical fit of the data to 
an infinite number of models we can just take a picture 
of the sky and from this obtain a picture of correlated 
maps. As already argued, the spectrum of the spots can 
be measured to judge if we are really looking at the topo- 
logical lensing of the horizon at the time of decoupling. 
The idea of combing the data for circles in the sky also 
shares the model independence feature and motivation. 

We apply our method to two small hyperbolic topolo- 
gies, the Weeks space and the Best space, named after 
their discoverers. The Best space is a compact hyper- 
bolic manifold obtained by identifying the twenty faces 
of a regular icosahedron |24|. The Weeks space Q has a 
more complicated fundamental domain with 18 faces and 
is of particular interest since it is currently the smallest 
compact hyperbolic space known. 




FIG. 21. The Dirichlet domain for the Weeks space. 



A. Strong patterns in a Weeks space 

It is advantageous to consider H 3 embedded as a 3D 
surface in a AD Minkowski spacetime. The universal 
cover H 3 is then a pseudosphere and the AD coordinates 
are restricted to the 3D surface with pseudoradius —1, 



-1 



(17) 



The isometries of H 3 can then be written as 4 x 4 matrices 
of the special Lorentz transformations. The coordinate 
transformation 



u 

Ui 
U 3 



coshr 

sinh r sin cos < 
sinh r sin sin <; 
sinh r cos 



(18) 



recovers the cosmologically familiar form of the 3D met- 
ric distance in comoving coordinates, 



ds 2 = dr 2 + sinh 2 r(d0 2 + sin 2 



(19) 



The geodesies take on a particularly simple form |26|] 
in the Minkowskian space. For simplicity, we tended to 
leave the Earth at the origin. For completeness, we men- 
tion how to move the Earth away from the center of the 
universe. If we align the Earth with the z-axis at position 
u e ~ ( u °i 0, 0, ul) we can Lorentz boost the Earth to the 
origin with the transformation 



A 



/ 7 -7/? \ 
10 
10 
7 / 



(20) 



and 7 = u° e and j/3 = u^. A photon observed in the 
direction n has v = — h and thus originated on the SLS 
at coordinates 



/ coshA^ 
nsinhA?7 



(21) 



This can easily be generalized to an arbitrary location 
with an arbitrary Lorentz boost. 
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The radius of the SLS, A77, depends on the value of Q t 
and the redshift of last scattering. In general, 



r\ = arccosh I 1 + 



2 - 20, 



(22) 



0,(1 + *) 



in units of the curvature radius. The volume of the SLS 
is 



The volume of the SLS grows exponentially with time. 
Therefore, many more copies of the fundamental domain 
can be contained within a surface of last scatter. The 
number of copies is a topological invariant, quite unlike 
in flat space. In flat space, one is free to set the volume of 
the manifold relative to the volume of the SLS arbitrarily. 
This is not possible for H 3 , as is ensured by the rigidity 
theorem [ p7| , which states that the volume in units of 
the curvature is fixed. A peculiar consequence is that the 
volume of a manifold is a topological invariant. In effect, 
if we measure topology we can refine our measure of the 
radius of the last scattering surface and thereby ft Q . In 
other words, we can use topology to measure curvature. 

In order to compactify the space, we consider the spe- 
cific example provided by the Weeks space. The fun- 
damental domain is a polyhedron with 18 faces and 26 
vertices shown in Fig. gl| which was taken from Snap- 
Pea, a census of compact hyperbolic manifolds j28| . The 
Weeks space is the smallest 3-manifold known with a vol- 
ume of ~ 0.94. With o = 0.3, and the redshift of last 
scattering taken to be z = 1100, so A77 = 2.328, we get 
Vsls = 150.64, and so there are roughly 150 copies of 
this universe within the SLS. With ft Q ~ 0.6, An ~ 1.5 
and there are only about 5 — 6 copies. 

There are 9 identification rules to glue these 18 faces in 
pairs. The 9 gi can be used to define a set of generators 
with many relations among them (such a set of genera- 
tors was also used by Fagundes Q ) . They are related to 
the words a, b of a much simpler presentation of the fun- 
damental group {a, b : ababa~ 1 bba~ 1 b, abab~ 1 aab~ 1 ab}. 
We extract from SnapPea these nine face-pairing 0(3, 1) 
matrices and their inverses out to 12 significant digits. It 
is necessary to use extremely precise matrices due to the 
chaotic flows. The sensitivity to initial conditions causes 
the image points to quickly be lifted off the pseudosphere 
if insufficient precision is used. 

The image points are denoted by in Minkowski co- 
ordinates: 




7r(sinh(2A77) - 2An) . 



(23) 



n 



and hi = ±1, .., ±9. Since we number the g^ from 1 — 9, 
let g Q be the identity so that yft = u /fI . We want to find 
the closest separation between and any of the image 



points y^. The gki---9k n form words of length n where n 
is the highest order of the farthest neighbor. 

Another way to count copies is to count periodic 
geodesies. The identification rules correspond to the min- 
imal closed loop geodesies. The number of geodesies of 
length L goes as 

N ^ ~ HL 

which is a result from chaos theory, with the Kolmogorov- 
Sinai entropy determined by the volume scale, h ~ 
F -1 / 3 . If we assume crudely that L = nV 1 ^ 3 so that 
n is the order of the neighbor, that is, the length of the 
word y^, then 

N(n) ~ e n /n 

So, for Vt = 0.3, the SLS encompasses roughly 150 copies 
of the fundamental domain and includes copies which 
are between 6 and 7 words away. For Q Q = 0.6, the 
farthest neighbors are about 3 words from the origin. 
The 9 generators and their inverses combine to form 18 n 
words of length n but as a result of the many relations 
among the all but N(n) ~ e n /n of these are repeats. 
Scanning the 18 n possible images N P i X x N v i x times is 
a huge numerical demand. To manage the task we use 
the following numerical algorithm, which is analogous the 
scheme used in the flat cases: 

(1) Move the point u on the SLS in towards the origin 
two steps to the image points y% 2kl = gk 2 9k 1 ^ tJ '- Out of 
these 18 2 points, select the three nearest the origin. 

(2) Move these 3 images two more steps and keep 
the one image nearest the origin so that it is essen- 
tially within the fundamental domain. Depending on 
how many neighbors fit within the SLS, step (2) can be 
repeated. 

(3) Repeat the same steps to the point u' being com- 
pared. Both points should now be located within the 
fundamental domain. 

(4) Keeping one image fixed within the fundamental 
domain, move the other within two neighbors until the 
distance between the two images is minimized. 

(5) The geodesic distance on the pseudosphere is 

d(y^y^) = arccosh(^V° ~ V ' y) • (24) 

We use this in our approximation 

C M (n,h') « C u {d^{x{n),x'{n'))), (25) 

where we use the angular correlations obtained from the 
CMBFAST code, with Vt h = 0.05 and the remainder of 
the subcritical density is assumed to be made up of dark 
matter. 

It is important to note that in a negatively- curved cos- 
mos there is, in addition to the Sachs- Wolfe effect on 
the surface of last scattering, an integrated Sachs- Wolfe 
(ISW) contribution to the perturbations. The ISW ef- 
fect is due to the decay of the gravitational potential as 




FIG. 24. 



the photons transit the space. Although we have not yet 
fully included the ISW effect, only the fluctuations on 
the largest scales should be affected whereas spots probe 
small-scale physics. On those scales where the ISW ef- 
fect contributes the geodesies are deviating sufficiently so 
that photons that originated in the same vicinity quickly 
take different paths with different decaying gravitational 
potentials. Correlations will not therefore be enhanced. 
One might fear that topological correlations could be 
erased by the different histories of the two initially adja- 
cent photon trajectories. Since the ISW effect does not 
effect the Doppler peaks in an infinite cosmos, it should 
leave the universe's spots unmarred. The next phase of 
investigation should include the ISW effect. Our maps 
correspond to data for which the lowest multipoles have 
been cleaned off. 

The first correlation we compute compares a random 
point with the rest of the surface of last scattering, 
C(u^, n), as shown in Fig. ^2|. This is unlike any correla- 
tion we have considered so far, but it can clearly be quite 
successful at uncovering geometric properties. As in the 
flat universes, the spots are likely spread out since we 
have not accounted for the discretization of the harmon- 
ics of the finite box which causes big dips in power at large 
modes, especially for such a small space. Again, the in- 
clusion of a huge number of terms would be needed to in- 
corporate this effect. We also compute C{fi,g%n/\\g%n\\) 
in Fig. H|. This transformation combines a boost along 
the hyperbolic surface with a rotation. 

B. The Best patterns 

Best built three compact hyperbolic spaces by identi- 
fying the faces of an icosahedron possessing twenty trian- 
gular faces, as shown in Fig. |24|. All of these spaces have 
volume V = 4.6860342738. The smallest geodesic ball 
which can encompass this space has a radius of radius 
r ~ 1.38 while the largest sphere which can be inscribed 
within the icosahedron has radius r ~ 0.87. For a rough 
estimate of the number of copies we count N = Vsls/V. 
With Q G = 0.3, there are on the order of 32 copies of this 
Best space within the SLS. 
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FIG. 25. The correlation of one point on the surface of last 
scatter with the rest of the sphere. The point is near the 
origin of one of the clones of the fundamental domain. 




FIG. 26. The sphere correlated by the go transformation. 

The faces are identified with 10 generators, the matrix 
representations for which can again be found in Snap- 
Pea. The abelianized homology group is Z/35. Fagundes 
also studied this space and, in particular, the occurrence 
of periodic quasar images [Q. The other three non- 
isomorphic Best spaces with icosahedra as fundamental 
domains have different fundamental groups and homol- 
ogy groups of Z/29 and Z/2 x Z/2. To construct our 
maps, we follow the same procedure as detailed for the 
Weeks space. The cosmic soccer ball in Fig. ^5] is the 
correlation of one point on the SLS with the rest of the 
sphere. The point happens to be very near the origin of 
one of the copies of the icosahedron. The patterns in the 
plot reflect the extreme symmetry of the fundamental do- 
main and also hint at the fractal nature of the geodesies. 
Notice the five-pointed star surrounding the tetrahedron. 
In the triangular corners of the five-pointed star there ap- 
pear to be six-pointed stars surrounding hexagons. We 
also show an example of the correlation C(n, geh/\ \gefi\ |) 
in Fig. ^6|. 

As data from the planned satellite missions becomes 
available, the CMB can be scanned for any hidden ge- 
ometric features. In the meantime, these correlated 
spheres show the huge potential for a pattern-oriented 
search of topological lensing. 

Our real-space approximation allows one to calculate 
temperature correlations while avoiding the analytically 
intractable eigenvalue problem on compact hyperbolic 3- 



spaces. Once a candidate compact universe has been es- 
tablished, such templates are useful for deeper statistical 
studies ||. However, an observational search for topol- 
ogy through large-angle temperature correlations can be 
made without a particular template being presupposed. 
This model independence of a pattern-driven approach is 
particularly important because it allows us to take pic- 
tures even if the space does not have constant curvature, 
if observations depend strongly on the location of the ob- 
server, or if chaotic mixing leads to unusual primordial 
spectra. Since we cannot predict topology within existing 
theories, we need this flexibility. 

While general relativity predicts the evolution of cur- 
vature, it does not specify the topology of space or of 
spacetime. Only a theory beyond Einstein's will be able 
to fully specify the geometry of the universe. Supergrav- 
ity theories necessarily acknowledge the importance of 
topology and compact hyperbolic cosmologies have even 
recently been studied as a consequence of string theories 
[ |3l| . However the universe was born, it was endowed with 
some topology and will have a place in our cosmic zoo. 
Astronomical studies of the topology of the universe may 
provide the most important insights into those aspects of 
the fundamental laws of Nature that dictated the global 
character of the space in which we live. 

Geometry effects not just the large-scale universe, but 
also the animals which inhabit it. As progeny of the 
universe, animals inherit certain cosmic blueprints. Per- 
haps it is only fitting that many analogues to the cosmic 
patterns could be found here on Earth on the backs of 
insects, in animal markings, even in human made mon- 
uments. If we could create a zoo of universes, each with 
a different topology, we might replicate all the animal 
markings from zebra stripes to leopard spots. 
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